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We use the proximity effect to generate effective topological superconductors by placing metals 
with strong spin-orbit coupling in contact with a superconductor, aiming to produce Majorana 
zero modes useful for topologically-protected quantum computation. In recent experiments, several 
quintuple layers of Bi2Se3 were epitaxially grown on the high-Tc material Bi2Sr2CaCu208-i-5, and 
conflicting experimental results were reported. We use the standard mean-field approach to study 
this heterostructure and find it is unlikely to have a large proximity-induced superconducting gap. 
Despite the seemingly correct temperature dependence, the s-wave gap claimed to be observed may 
not be purely superconducting in origin. Future work on the proximity-induced bulk superconduct¬ 
ing gap and the interfacial bandstructure should shed light on this issue. 

PACS numbers: 74.20.Rp,74.45.+c,74.78.Fk 


I. INTRODUCTION 

The proximity effect in superconductor-metal inter¬ 
faces, or the leakage of Cooper pairs from the super¬ 
conductor into the metal, has been studied since its 
experimental discovery in 1960 ^. The magnitude of 
the induced superconducting gap generally depends on 
the quality of the interface, external fields, and certain 
material-dependent properties, in particular the Fermi 
surfaces and energy spectra of both the superconductor 
and metafSTi^. The symmetry of the induced supercon¬ 
ducting gap is typically identical to that of the supercon¬ 
ductor unless there exists symmetry-breaking perturba¬ 
tions, such as the lattice-symmetry mismatch, ferromag¬ 
netism in the metal^i^, or the presence of large spin-orbit 
interactions^. 

Recently, the application of the proximity effect to 
the surface states of a topological insulator (TI), with 
the goal to produce the topological superconductiv¬ 
ity, has stimulated intense theoretical and experimen¬ 
tal effort. In addition to s-wave superconductors 2 r— , 
d-wave superconductors also show interesting proxim¬ 
ity effect on the surface states of a Tl^idirJ^. Wang 
et alii^ observe a surprisingly large s-wave-like pair¬ 
ing gap on the top surface of a seven-quintuplc-layer- 
thick Bi2Se3 film on Bi 2 Sr 2 CaCu 208+5 (BSCCO) by us¬ 
ing angle-resolved photoemission spectroscopy(ARPES). 
By applying a suitable incident photon energy (approx¬ 
imately 50 eV), at which the contribution from the sur¬ 
face states dominates that from the bulk states, a gap 
of up to 15 meV on the outermost surface states of the 
Bi2Se3/BSCCO system is detected. The change of the 
pairing symmetry from d-wave in BSCCO to s-wave-like 
on the surface states from the proximity effect is claimed 
to reflect a nontrivial coupling in this systemic. However, 
two separate experimental group a^^d^ with similar setups 
do not reproduce the proximity-induced pairing gap on 
the Bi2Se3/BSCCO. Earlier tunneling measurements^^ 
suggest the proximity induced gap at the interface is d- 
wave. These conflicting observations prompt us to ad¬ 


dress the issue theoretically and ask whether such a large 
proximity-induced s-wave gap is possible. 

In this paper, we use the tunneling model pioneered 
by McMillanii to study the superconducting proximity 
effect between Bi2Se3 and BSCCO. We assume the inter¬ 
face between the superconductor and the metal is clean 
or with only weak random impurities, and the super¬ 
conducting gap in BSCCO is identicaU^ for all coupling 
strengths between BSCCO and Bi2Se3. We model the su¬ 
perconductivity within BSCCO via mean-field, and eval¬ 
uate within Bi2Se3 the induced pairing amplitude, de¬ 
fined as the sum of the expectation value of Cooper pairs 
in momentum space, for a range of symmetries. Our re¬ 
sults show that the proximity-induced superconducting 
gap is small and mainly d-wave, and the unambiguous 
way to find the condensates is through the evaluation of 
pairing amplitudes. The smallness of s-wave and p ± ip 
pairing amplitudes supports the absence of a large in¬ 
duced s-wave-like superconducting gap. 

Based on our results, we suggest the discrepancies be¬ 
tween the experimentsi^^— are due to different interface 
coupling strengths between Bi2Se3 and BSCCO across 
the various samples. The observed gap in ReflU may 
not purely be due to superconductivity, but also to a 
mass gap generated by bands crossings caused by large 
interfacial tunnelings as compared with the samples used 
in Refill and m The growth temperature in Refill is 
certainly much higher than that in RefUfil suggesting the 
interface tunneling strength is larger in sample of Refill. 
Further analysis on the symmetry of the bulk sup ercon- 
ducting gap of Bi2Se3 in the samples of Refill should 
reveal its d-wave nature, as the bands overlap between 
the bulk bands of Bi2Se3 and BSCCO should be smaller 
than the overlap between the surface states. Another 
way to verify our claim is to measure the band structure 
of Bi2Se3 at the interface, which should show a d-wave 
single particle gap structure with lobes pointing in the 
nodal direction when the system temperature is higher 
than Tc, the superconducting transition temperature of 
BSCCO. 
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We note that earlier work by Z. X. Li et. alj^^ has 
also addressed the same issue with a two-bands model of 
the Bi2Se3. In their work there are two quantum well 
states other than the surface state present at the Fermi 
surface for the 7 QL-thick Bi2Se3 . They also introduce 
small (around 100 meV) repulsive and attractive inter¬ 
actions within the Bi2Se3. At the top layer (away form 
the interface with BSCCO) they find s-wave dominant 
pairing due to the suppression of d-wave pairing near the 
r point and claim that disorder could further enhance 
s-wave pairing. The proximity induced superconducting 
gap magnitude in their paper— is also around 1 meV with 
similar tunneling str eng th. Their results support the re¬ 
ports in RefEi and [isl of the absence of a superconud- 
cting gap despite its s-wave gap structure. Furthermore, 
the dominant s-wave pairing gap in ReflTol is the result of 
multilayer structure of Bi2 803 (through enforcing s-wave 
interactions in Bi2Se3 and disorder, which prefer s-wave) 
and choice of Fermi surface around F point. This mech¬ 
anism is different from our work where the formation of 
s-wave-like gap is the combined effects of superconduct¬ 
ing proximity effect and hybridization of bands between 
Bi2Se3 and BSCCO. 

We organize this paper as follows: In SectionllTl we in¬ 
troduce our model Hamiltonians for the BSCCO/Bi2Se3 
system. Assuming good contact surface between the two, 
we compute various paring amplitudes and proximity- 
induced gaps in the surface state of Bi2Se3 at different 
tunneling strengths. In Section lHll we discuss other pos¬ 
sible factors which could enhance the proximity induced 
superconductivity. In Section llVl we find the combined 
effect of strong hybridization of bands and superconduct¬ 
ing proximity effect gives rise to s-wave-like gap struc¬ 
ture. In SectionlV]we summarize our results and suggest 
possible experimental approaches to verify our claim. 


II. PROXIMITY EFFECT: BSCCO AND BLSeg 


We start with the model Hamiltonian describing the 
Bi2Se3, a type of three dimensional topological insu- 
lator(TI), grown on top of the BSCCO. We assume 
only single metallic layer in the Bi2Se3 in our com¬ 
putation, which can be viewed as one of the two di¬ 
mensional surfaces of the three dimensional topologi¬ 
cal insulator Bi2Se3 in contact with BSCCO if the bulk 
band were insulating. The motivations for forming this 
heterostructure^iiir— is to produce effective p zL ip su¬ 
perconductivity, a form of topological superconductivity, 
on its interface. The edge modes or vertex states of the 
topological superconductor host the Majorana modest, 
the simplest anyon which could be useful for quantum 
computations. 

To describe the superconducting proximity effect be¬ 
tween surfaces of BSCCO and Bi2Se3, we use the follow¬ 


ing single band bilayer model Hamiltonian: 

H = HbSCCO + + Ht 

Hbscco = ^ dijEd{k) - 


k ,(7 




k,k' 


-ffBiaSea — ^ ^ ct^(Ac(fc)a,3 

k,oi,^ 

Ht='^(trcl dt -I -ct ). 


( 1 ) 


Here Hbscco^ HBhSes^ and Ht describe the Hamil¬ 
tonian for two dimensional surface of BSCCO, surface 
state of Bi2Se3, and the single particle tunneling terms 
between the two surfaces. Her, a, /3 denote spin in¬ 
dices, ^ and ^ the electron annihilation operators 
of BSCCO and Bi2Se3 in momentum space representa¬ 
tion, and k = {kx,ky) is the linear momentum of the 
two dimensional surface. lajS is the identity matrix, and 
Ec{k)ctp is a 2 X 2 matrix characterizing the spin orbit in¬ 
teraction of surface state of TI. Attractive V(k, k') with 
d wave symmetry is assumed to give d-wave supercon¬ 
ductivity of BSCCO under the BCS mean field approx¬ 
imation. The tunneling amplitude tk is proportional to 
the wavefunctions overlap between that of BSCCO and 
Bi2Se3. We have assumed few or weak random nonmag¬ 
netic impurities such that spin and momentum are con¬ 
served in the tunneling term, and constant tunneling am¬ 
plitude nearby the relevant Fermi level {tk = t). 


E/h 




(b) FS with warping 


FIG. 1 . (a) Energy dispersion of surface state of Bi2Se3 de¬ 
scribed in Eq. Q. The relevant energy scale in the main text is 
colored in shaded blue region, (b) We plot the Fermi surfaces 
of Bi2Se3 including the warping cubic terms at p.2/h ~ 1.3 
(red thin line) and p.2/h ^ 2.6 (blue thick dashed line), and 
BSCCO (green thick line). The ^2/d ~ 2.6 is just for illus¬ 
trating the hexagonal symmetry of the Fermi surface which 
is difficult to see with the Fermi level {fk2/h ~ 1 . 3 ) set by the 
experiments. 

The explicit form of Ed{k) and Ec{k)ai3 are obtained 
from the tight binding models, with parameters specified 
from first principle studies for Bi2Se3^‘2i and BSCC02^. 
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The single band dispersion for BSCCO is: 

Ed{k) = Ed{k) - ^d = -h (cos(fcxai) + cos(fcyai)) 

+ h' cos(fca;ai) cos(fcyai) — h (cos(2fca;ai) + cos(2fcyai)) 
— h (cos(2fca;ai) cos(fcyai) + cos{kxai) cos{ 2 kyai)) 

+ h (cos(2fc^ai) cos(2feyai)) — (2) 


Here h = 0 . 2975 eV, h' = 0 . 1636 eV, h" = 0 . 0259 eV, 
h!" = 0 . 0558 eV, h!'" = O.OSlOeV, Hd = - 0 . 1305 eV, 
oi ~ 3.82A. For the Bi2Se3 crystal, to account for the 
Cq crystal symmetry and different orientations on the 
proximity effect seen in the experiments, we include the 
warping terms^S in our model Hamiltonian. From Refl^. 
the low energy dispersion of Bi2Se3 close to the Brillouin 
zone center is given by 

Ec{k) = Ec{k) - fic 

= Po +Piikl + ky) +p2{cFxky - cFykx) 

+ P?,{k\ E k^_)az - Pc ( 3 ) 

where po = 0.035 eV, pi = 1.38556 eVA^, p2 = 0.795 
eVA, p3 = 0.3535 eVA^, k± = kx E iky, and ax, ay, 
and az are the 2 x 2 Pauli spin matrices. The p2 term 
is the Rashba spin orbit term. For our calculations, we 
need to extend this low energy Hamiltonian to the whole 
Brillouin zone in momentum space and keep the single 
Dirac-cone-like structure. The quadratic k terms have 
undesired behavior in the large-A: region and does not 
give the consistent band structures compared with exper¬ 
iments at large momenta. To include the warping terms 
and avoid the issues of inconsistency at large k, we keep 
the linear term unchanged, remove the quadratic terms, 
and replace kxa2 with sin(fca;a2) and kya2 with sin(A:ya2) 
for the third order warping terms. Here 02 = 4 . 138 A is 
the lattice constant of Bi2Se3. The low energy dispersion 
in Eq.(I 3 |) is then modified as: 

Ec{k) = Po + P 2 {(yxky - aykx) ( 4 ) 

-I- ^ {2sm^{kxa2) - 6 sm{kxa2) sin^ {kya2)) az - Pc 
0-2 


The chemical potential p^ exhibiting surface states of 
Bi2Se3 is estimated to be around 0 . 4 eV for the relevant 
experimenti^. Noting that Eq.Q is not of tight bind¬ 
ing form but from an effective low energy Hamiltonian 
of surface state of Bi2Se3 with parameters fixed by first 
principle studie a^°i^^ . The first order derivative of Edk) 
is discontinuous at zone boundary {kxa 2 , kya 2 being ± 7 r) 
after imposing periodic boundary condition. For the low 
energy limit (energy scale less than ^2^2 3 . 2 eE) dis¬ 

cussed in this paper this artifact does not enter into our 
equations or modify our results. The region of relevant 


energy scale in our discussion here is plotted in Fig 1 (a) 


For the numerical computations shown below we take 
nearest neighbor hopping h of BSCCO in Eq.Q as the 
energy unit. We use kxOi = k^ and kyUi = ky as dimen¬ 
sionless momentum space parameters although there is a 


minor difference in their lattice size (02/01 ~ 1 . 08 ). In 
this unit the d-wave superconducting gap is defined as 

k' 

= A*(cOs(fca:) — COs(fcy)) ( 5 ) 

Here the gap magnitude |A| ~ 40 meV (or |A|//i ~ 0 . 13 ) 
is assumed to be fixed (termed the ’’bulk limit” below) 
by viewing the superconductivity of BSCCO not influ¬ 
enced by the contact of Bi2Se3. The expectation value 
is taken with the ground state of the whole mean field 
Hamiltonian H. In principle we may also compute this 
gap magnitude self consistently, as shown for the model 
calculations in Appendices ETICl by fixing the interac¬ 
tion strength V{k,k'). The obtained proximity induced 
superconducting gaps or pairing amplitudes are smaller 
compared with those obtained with the bulk limit. To 
obtain an upper bound on the proximity induced gap on 
Bi2Se3 we stick with this bulk limit in the main text of 
the paper. 

We identify the proximity induced superconductivity 
in Bi2Se3 by computing two quantities: the pairing am¬ 
plitude and the quasiparticle energy gap. The pairing 
amplitude Ai is defined as: 


k' 


( 6 ) 


with fi being the symmetry factor and i denoting s, p, 
pEip, or d wave symmetry (a and /3 are spin indices, fs = 
1 , fp = sm{kx) or sin(fcy), /p±jp = sin(A:a;) ± ism{ky), 
and fd = cos{kx) — cos{ky). /' has similar definition 

as fi with k replaced by k'. This fif- geometric fac¬ 
tor comes from the angular expansions of V(k — k').) 
and Ai is the corresponding magnitude. Even (spatial) 
wave symmetry has odd spin angular momentum and 
vice versa, as required by the Pauli exclusion princi¬ 
ple. This dimensionless pairing amplitude directly re¬ 
flects the amount of Cooper pairs formed at Bi2Se3, but 
it is not directly probed by tunneling measurement^^ nor 
spectroscopy like ARPESi^^— . The physical quantity 
probed in these experiments is the quasiparticle energy 
gap. We identify the energy gap due to proximity ef¬ 
fect in two approaches: One is from the density of state 
(DOS) of Bi2Se3, and another is obtained from the nu¬ 
merical results of energy bands difference in the diago¬ 
nalized bases of total Hamiltonian H. The DOS calcu¬ 
lation is done with Hartree approximation (ignoring the 
exchange term), as shown in Appendix IA 11 and perform 
the momentum integral on the imaginary part of electron 
Green’s function of Bi2Se3. From the shape of the DOS 
we may identify the gap symmetry and magnitude, but 
will not have momentum space resolution. The numer¬ 
ically obtained energy bands difference with momentum 
dependence complements this. 

From Eq.Q we find the desired surface state band 
structure of Bi2Se3, showing clear hexagonal Fermi sur¬ 
face at larger chemical potential as shown in Fig|l(b)[ 






4 


t 


Ap 

Ad 


0.01 

-7.1698 

-0.29080 

2.7588 

-2.9932i 

0.1 

-8.7897 

-4.9422 

209.67 

-47.220i 

0.67 

-364.21 

-93.741 

8713.9 

-1254.4i 

1.0 

-545.43 

12.993 

17739 

-1381.5i 

1.5 

-406.05 

768.48 

32251 

1150.21 

2.0 

1475.2 

2659.3 

41063 

6947.71 


TABLE I. Various possible pairing amplitudes at different 
tunneling strengths t. Gap magnitude of BSCCO is fixed at 
A//i ~ 0.13 with pairing amplitudes evaluated with 160 x 160 
k points. Ai is defined as Ai = Ai x 10“®. 


This hexagonal structure is not apparent at the exper¬ 
imental relevant Fermi level. However, with the break¬ 
ing of circular symmetry, the presence of Zeeman like Uz 
terni^^, and the crystal orientation of Bi2Se3 which is 45 ° 
or 15 ° differentia from that of BSCCO, we have generated 
s-wave, p-wave, and p ± ip-wave pairing in additional to 
d-wave pairing on the Bi2Se3 surface. The crystal orien¬ 
tation difference, as demonstrated in Fig |l(b)| is impor¬ 
tant for the existence of s-wavei^ pairing. There would 
be no s-wave pairing if rotated by 15 ° (or equivalent 45 ° 
with Cq symmetry in warping corrected Bi2Se3 surface 
state) due to the d-wave pairing from the BSCCO and 
the alignment of symmetry axes of Bi2Se3 in nodal di¬ 
rections. Same argument also applies to the lacking of 
s-wave pairing without the inclusion of warping terms. 
Aforementioned pairing amplitudes are denoted by Ag, 
Ap, Ap±ip, and Ad and their values with different tun¬ 
neling strengths t are listed in TablelH 
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FIG. 2. (a) DOS in Bi 2 Se 3 with t/h — 0.67 and |A|//i ~ 0.13. 
Blue line is the raw numerical data and purple line is the poly¬ 
nomial fit for easier identification of gap location. Frequency 
unit is meV. (b) Angular dependence of gap in momentum 
space obtained from the bandstructure in Bi 2 Se 3 /BSGCO sys¬ 
tem with warping terms, t/h = 0.67 with angular separation 
of 5°. 

In Fig|^ we plot the DOS of Bi2Se3 surface state at 
tunneling strength t/h = 0.67 with superconducting gap 
A//i ~ 0.13 in BSCCO. The proximity induced gap is 
mainly d-wave-like as shown in TableHl with d-wave gap 
magnitude around O.OOI6/1 or 0.5 meV (about half from 
DOS plot due to the d wave nature). This magnitude is 
slightly larger— than the gap obtained from the difference 


of particle hole bands of Bi2Se3 nearby the Fermi level 
shown in Fig | 2 (b)[ which is around 0.3 meV evaluated 
with 160 X 160 k-points. The s-wave-like component read 
out from Fig | 2 (b)| is about 0.08 meV, which gives roughly 
consistent ratio as (Ag + Ap±ip)/Ad shown in TableUDthe 
p ±ip also gives s-wave-like gap). 

In the ARPES experiment of Refill seven quintuple 
layers (QL) of Bi2Se3 are used, and the measurement 
is done at the top plane which is not in contact with 
BSCCO. The model calculations done above is just for 
a single (surface) layer of Bi2Se3 in direct contact with 
BSCCO, assuming the thickness of Bi2Se3 is sufficient 
large that the Dirac modes on the two opposite surfaces 
are present. As we have assumed a very smooth con¬ 
tact between the two materials, the gap magnitude ob¬ 
tained, via DOS or rediagonalized band structures, in 
our calculations (corresponding to the proximity induced 
gap at the surface in contact with BSCCO) should be 
viewed as the maximum value. Our result shows that 
the gap is small (< 1 meV at t/h = 0 . 67 ) and mainly d- 
wave like even with the inclusion of warping and correct 
crystal orientation, which tends to support the absence 
of the proximity effec t in the surface states of the 
Bi2Se3/BSCCO system. Below we discuss other factors 
not included in the aforementioned discussions, and see 
how they modify our results. 


III. OTHER RELEVANT FACTORS 

We have assumed smooth interface and computed only 
the proximity induced gap at the surface of Bi2Se3 in di¬ 
rect contact with BSCCO. In the actual experimentsi^^— 
the ARPES were done on the opposite surface of a few 
quintuple layers thick Bi2Se3. From the ARPES data the 
Fermi level of a few QL thick Bi2Se3 actually crosses the 
bulk band and the bulk is not insulating. Furthermore, 
in the growth process, there are always some impurities 
generated at the bulk or surface of Bi2Se3. Last but not 
least, the lattice size of BSCCO is S.SA and that of Bi2Se3 
is 4 . 138 A. The latter is taken to be the same as that of 
BSCCO in our mean field calculations. We discuss how 
these factors influence our results below. 

The simplest way to consider the proximity effect on 
the other surface of Bi2Se3 is to assume the two sur¬ 
faces are tunnel-coupled. In this case the model becomes 
trilayer rather than bilayer one. For moderate tunneling 
strength the multilayer structure always has smaller pair¬ 
ing amplitude compared with bilayer system (as shown 
in the generic model discussion in AppendixlB]). Thus 
the simplest way to include this finite thickness effect of 
Bi2Se3 will not be able to account for the 15 meV mea¬ 
sured in ReflU 

We may include the bulk band by treating it as a two 
bands metal on the metallic surface in contact with the 
BSCCO. This is because the Fermi surface of the Bi2Se3 
also cut through the bulk band, and in some literature^ it 
is called a topological metal. In principle the bulk bands 
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(a) DOS (b) Comparison 

FIG. 3. (a) The DOS of Bi 2 Se 3 at tjh — 2 with supercon¬ 
ducting gap in BSCCO A//i ~ 0.13. All other parameters 
are the same as those in the main text. Frequency in unit of 
meV. (b) Gap structure (red dots) at t//i = 2 compared with 
the experimental data (black square) shown in Refill. 


discussion in AppendixlC]) . However, given the small lat¬ 
tice mismatch (lattice size ratio 02/01 ~ 4 . 1 / 3.8 ~ 1 . 08 ) 
the major enhancement of proximity induced gap should 
come from the enlargement of Fermi surface on the 
Bi2Se3 surface. Based on our calculations, even rais¬ 
ing to twice of the original chemical potential of Bi2Se3 
does not significantly enhance the proximity induced gap 
magnitude (around twice and also mainly d-wave-like). 
Judging from the discrepancy of shape and gap magni¬ 
tude with ARPES data shown in Refill, and the fact 
that the estimates we have here should be the theoreti¬ 
cal upper bound for a clean interface, we tend to support 
the experimental results in RefsIU andfisl 


share the tunneled Cooper pairs from the superconductor 
with the surface band, therefore at the same tunneling 
strength the pairing gap at the surface decreases com¬ 
pared with that without the inclusion of bulk band in a 
self consistent calculation^^. If the two bands are cou¬ 
pled via repulsive interactions, the pairing gaps at differ¬ 
ent bands could be different as suggested to be measured 
experimentally in Refill, where different frequencies are 
used to distinguish the measurement of bulk from that 
of surface. As both bands share the superconducting gap 
from BSCCO, the inclusion of bulk gap tends to lower 
the proximity induced gap of the surface band, rendering 
it even smaller than 0.6 meV at tjh = 0 . 67 . 

The impurities in the bulk and surfaces of Bi2Se3 could 
in general suppress the nodal pairings, and make s-wave 
pairing dominant for sufficient thick sample^^. The draw¬ 
back is the scattering from the impurities also tend to 
diminish the pairing gap, making it unlikely to achieve 
the large (> 10 meV) proximity induced s-wave-like gap. 
Another way to enhance this s-wave component is to in¬ 
troduce an isotropic, attractive interaction as in Reffiol. 
However, in our zero temperature calculations, including 
this attractive interaction will make the Bi2Se3 super¬ 
conducting without coupling to BSCCO. The undoped 
Bi2Se3 has not been seen to be superconducting in the 
experiments^, and therefore we do not include this fac¬ 
tor in our discussion. It is possible that interfacial s 
wave superconductivity could be formed at the inter¬ 
face of BSCCO and Bi2Se3, and we verify this interfa¬ 
cial attractive interaction does enhance the s-wave pair¬ 
ing amplitude in Bi2Se3 layer but the enhancement is 
small (around 25 %) with moderate interaction strength 
(twice the attractive interaction strength of BSCCO at 
tjh = 1 . 3 ). 

There are two other factors which could possibly en¬ 
hance the pairing amplitudes in this system. One is 
the lattice mismatch issue, and the other is the possi¬ 
bility of forming enlarged Fermi surface in Bi2Se3 due to 
its growth on BSCCO (the gating issue). For different 
Fermi surface structure the lattice mismatch could usu¬ 
ally lead to more overlaps in bands when taking the lat¬ 
tice mismatch into considerations (see the generic model 


IV. POSSIBLE S-WAVE GAP 

How do we explain the experimental findings in Refill 
It is suggestedi^ that the grounding issues could lead 
to incorrect readings of superconducting gaps, but this 
factor might not reflect the temperature dependence of 
the gap in a way correlated with the BSCCO observed 
in Refflg We suspect the observed s-wave-like gap is 
the result of strong hybridizations between the bands of 
Bi2Se3 and BSCCO due to the large tunneling ampli¬ 
tude. In this case the observed gap is not purely due to 
proximity induced superconductivity, but also the mass 
gap generated by single particle tunneling. 

We find for tjh > 2 the single particle mass gap along 
all directions (both nodal and antinodal ones) is gen¬ 
erated even if we turned off the superconductivity in 
BSCCO, and we use this as the upper bound^l of tun¬ 
neling amplitude as no gap opening along antinodal lines 
is found at T > Tc in Reflla. From Table U we find at 
t/h = 2 the proximity induced pairing amplitudes are 
still dominated by d -wave pairing, but the shape of the 
gap as shown in Figj 3 (b)| is similar to s-wave. The gap 


is around 13.4 meV at t/h = 2 , and no gap along the 
antinodal direction of BSCCO is seen when the super¬ 
conductivity is turned off. The lobes of this mass gap is 
pointing along the nodal directions of BSCCO, making 
it look like s-wave gap structure, as shown in Fig | 3 (b)[ 
when the superconductivity is turned on. 

The aforementioned s-wave-like gap opening is con¬ 
sistent with the temperature dependence observed in 
ReflU, i.e. no gap at T > Tc and s-wave-like gap at 
T < Tc- At t/h = 2 the onset of proximity induced super¬ 
conductivity is estimated, from the change of DOS, to be 
around 60 K which is roughly consistent with the datapi^. 
Factoring in the aforementioned lattice mismatch, en¬ 
larged Fermi surfaces, possible s-wave pairing interaction 
at the interface, could possibly give rise to the scale and 
features similar to the experimental results in Refill. We 
emphasize that the s-wave-like structure shown in Fig| 3 ] 
is the combined results of bands hybridization between 
Bi2Se3 and BSCCO and the d-wave superconductivity 
from BSCCO. The proximity induced pairing amplitude 


magnitude read out from the DOS of Bi2Se3 in Fig 3 (a( 
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in Bi 2 Se 3 is still dominated by d-wave. This scenario 
is different from the dominant s-wave pairing discussed 
in Refill, in which the s-wave pairing gap is completely 
from the proximity induced superconductivity. 


V. CONCLUSION 

In this paper, we use the mean-field approach to com¬ 
pute the superconducting pairing amplitudes and gap 
magnitude in a superconducting-metal interface. We 
find the pairing amplitude is mainly d-wave except at 
small t {t/h < 0.01 with h the nearest neighbor hop¬ 
ping amplitude in BSCCO). With reasonably large t 
{2 > t/h > 0.1), the gap structure, obtained from our 
DOS and band calculations, is mainly d-wave with gap 
magnitude less than 1 meV. 

We suggest one possible reason why the large s wave 
like gap (~ 15 meV) is observed in ReflU while the 
other two similar setup a^^d^ show no significant (> 5 
meV) proximity induced gap. We argue the gap in 
BefllOl is not purely due to superconductivity, but rather 
contains a large component due to a large tunneling 
amplitude(t//i > 2). The pairing amplitude evaluated 
at large t are still dominated by the d-wave channel, and 
its mixing with the single-particle gap makes it more s- 
wave-like, as shown in Fig |3(b)[ Thus, we propose that 
sample-to-sample variations in interfacial coupling be¬ 
tween Bi 2 Se 3 and BSCCO explain the discrepancy be¬ 
tween the observed gaps, and also that the induced su¬ 
perconductivity is much smaller than would be naively 
concluded from the observations of Refill. 

A piece of supporting evidence for this claim, as also 
pointed out in R efill and 111, is the lack of superconduct¬ 
ing coherence peak in the observed ARPES spectrum in 
Ref. Ill. Furthermore, it is claimedi^ that the bulk con¬ 
ducting bands of Bi 2 Se 3 are also superconducting, but 
with a much smaller gap. As the bulk bands are fur¬ 
ther away from the quasiparticle bands in BSCCO and 
their effective tunneling amplitudes are smaller due to 
smaller wave function overlaps, it is possible that the gap 
induced in the bulk bands is dominated by the d-wave 
superconducting gap and the measured gap might be 
mainly from proximity induced superconductivity. More 
detailed measurements of the bulk spectrum should shed 
light on this issue. This case also serves as a caution¬ 
ary example: measuring the temperature dependence of 
a gap at an interface is not sufficient to show that the 
gap is superconducting in origin. 
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Appendix A: Perturbative results for proximity 
induced pairing amplitude 


We derive the analytical, perturbative results of su¬ 
perconducting proximity effect for some generic bilayer 
and trilayer systems using Green’s functions formalism. 
The perturbative term is the attractive interaction in the 
superconducting layer, and we further assume the tunnel¬ 
ing amplitude is small to simplify our results. The form 
of self energy in the Green’s function is assumed to take 
the Hartree form (neglecting exchange interactions) with 
mean field approximations. The real part of self energy 
just shifts or renormalizes the chemical potential. Under 
these approximations the self consistent equation of the 
superconducting layer is the same as the uncoupled one, 
and the leading term of pairing amplitude in the first 
metallic layer is proportional to The generic model 
Hamiltonian we considered in the Appendix is: 


k,l,a,(7' 


+ ^ 






(Al) 


Here l,i,j range from 1,2,3,..., A with 1 denote elec¬ 
tron/hole operators on the superconducting layer and 2, 
3 (say A = 3) denote two bands in one metallic layer 
or two metallic layers. In this model Hamiltonian we 
consider the usual form of density density interactions. 
Below we show explicitly the calculations for three cases: 
Bilayer, trilayer, and metallic layer with spin orbit inter¬ 
actions. 


1. Two layers with single band 

Here we perform perturbative calculations via the 
Green’s function formalism. Using the four component 
Nambu spinor bases = (ct^^ 
the Green’s function is given by 

G(fc,r) = -(T,^(fc,T)«'t(fc,o)) (A2) 

For Vi{q) = 0 the retarded Green’s function Go{k,uj) 
takes the form: 
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Goik,uj) 


^ uj — ei(fc) + i 


IT] 


V 


0 

t 

0 


0 

w + ei(fc) + ir] 
0 

-t 


t 

0 

uj — e2(fc) + iv 
0 


0 V' 

—t 

0 

w + e2(fc) + irj J 


(A3) 


with e;(fc) = e;(fc) — /i;. Turning on Vi(^ 7 ^ 0 gives self 
energy correction on the bare Green’s function Go{k,uj) 
from perturbations in Vi{q). We take q = k + k' as in 
the BCS theory and denote Vi{q) = Vi{k, k') = Vi in the 
following. The perturbation in V) gives: 

G(fc, w)-i = Go(fc, w)-i - S(fc, uj). (A4) 

Based on first order perturbation in Vi{q), the Hartree 
corrected^^ l^{k,uj) ~ Ejy takes the following form 


^ Ei(fc,w) Si{k,u!) 0 

-Ei(-fc,-w) 0 

0 0 T,2{k,u}) 

V 0 0 


0 \ 

0 


S2ik,uj) 
-E2(-fc,-w) / 


with p{k,uj) = —^Q(G(k,ca)), fi{uj) = 9{pi — uj) the 
Fermi-Dirac distribution at zero temperature, and 


poo 

dujpii{9,u})fi{u}) 

k' 

^ Vn 

E2(A;,w)~^— / dujp33{k\u})f2{uj) 

^ J —00 
k' 

/ oo 

dujpi 2 {k',uj)fi{uj) 

-CX) 

k' 

/ oo 

dujp3i{k' ,uj)f2{uj) 

'OO 


(AS) 

(A 6 ) 

(A7) 

(AS) 


Ea. (|ASp -Ea. (|A8p is named Hartree approximation, and 
under this approximation the and Si terms have mo¬ 
mentum but no frequency dependence. For a given set of 
Vi, the Ei and Si can be obtained via iterations. We take 
El < 0 and E 2 = 0 to simulate the case of first layer be¬ 
ing superconductor and second layer being normal metal 
before placed together. Note that 


= I duJ^^iGil uj)i 2 ) (A9) 

duj-^^iG{k,uj) 34 ) (AlO) 


with G(fc, uj)ij denote i, j component of the 4x4 Green’s 
function. The momentum sum, along with appropriate 
symmetry factor fi dehned in Eq. ([5]) , in above equations 
give pairing amplitude mentioned in the main text. The 
DOS of metallic layer p{uj) is computed by 


1 _ f dk 

p(w) = -Yl 1™^'^ / prA2‘^*^(^ + *^’^)j-(Aii) 

Using assumptions Ea. (IA5l) -Ea. (IA8l) and the conditions 
Ul 7 ^ 0, V 2 = 0 in Ea. (|A4ll we have 


G(fc,C.)l2 = 
G{k,Uj)3A = 


_ Si{k,uj){uj'^ - e 2 (fc)^)_ 

— e2(fc)^)((ei(fc) -I- Ei(fc,a;)^ -I- S'i(fc,a;)^ — oj'^) + 2t'^{uj'^ -|- e2{k){ei{k) -|- Ei(fc, w))) — U 

—t‘^Si{k,uj) 

{uj^ — e2(fc)^)((ei(fc) -I- Ei(fc, a;))2 -I- Si{k, uj)^ — u)^) + 2t^(uj^ + e2(k)(ei(k) + Ei(fc, w))) — U 


(A12) 

(A13) 


with UJ = UJ + ip for retarded Green’s function. For t = 0 

we have (ct ^ ^ ,) = 0 and 
' fc 2 t -k 2 i' 


(ct C^ - ) = 
' fclt -fclj.' 


-S3{k,Uj) 


2J(ei(fc) + Ei(fc,w))2 + ^i(fc,u;)2 


-Si{k,uj) 

2Ei{k) 


(A14) 


The single particle self energy correction terms E^ can 
be absorbed into the renormalized chemical potential pi, 


and we remove this term in the following expressions for 
simplicity. 


By combining Ea. (IA14p with Ea. (IA7|) . we get the self 
consistent gap equation in the BCS theory by identify¬ 
ing superconducting gap Ai{k) = Si{k,uj). For small 
but finite t we take the leading correction in the eval¬ 
uation of frequency integral for G{k,uj) 34 ^. Under this 
small t approximation the pole structure is the same as 

















in Ea. (|A14[) . and we have: 


fe2t -fc2i 2Ei{k){Ei{kY - e2{ky) 

Note that it is straightforward to obtain full analytic re¬ 
sults without this small t (tunneling strength) approxi¬ 
mation, but the results are less illuminating. We stick 
with the Hartree and small t approximations in this Ap¬ 
pendix 1^ to illustrate the main idea. 

From Ea. (|A15|) the rough estimate gives the second 
layer superconductivity gap magnitude IA 2 I: 


IA 2 I 


ksTc^ < 


e 

\4-El\ 


|Ai| 


t^ksTci 

\4-El\ 


(A16) 


Here we use IA 2 I < I Lfe ^i(42t^-fe2i^l’ 

as the average of Ei(k) and ei{k) in momentum space 
as a rough estimate. From this estimate, we know for 

I 


( 00 

- ei(A:) -b ir\ 

0 

tl2 


0 

00 -b ei(fc) -b IT] 

0 


tl2 

0 00 

- l2(k) 


0 

—tl2 

0 


G 3 

0 

^23 

V 

0 

—tl3 

0 


e 2 (fc) very different from ti{k) (’’mismatched” Fermi sur¬ 
face) the denominator in Ea. (IA16l) increases, giving rise 
to smaller gap magnitude as expected. The gap symme¬ 
try of the second layer is determined by both the first 
layer gap symmetry and the Fermi surfaces from both 
layers. We do not use this estimate in the main text, 
but Eq. (IA16I) gives naive intuitions why the mismatched 
Fermi surfaces give smaller proximity induced pairing 
amplitude. 


2. Two bands or three layers 


For two metallic bands or two metallic layers we follow 
the same definition of generalized Green’s function and 
extend the 4x4 bases to 6 x 6 ones to accommodate 
this extra degree of freedom. For Vi{q) =0 we have the 
retarded Green’s function Go{k,uj) taking the form: 


0 

tl3 

0 

\ 

-J. 

— ^12 

0 

—tl3 



0 

t23 

0 


(A17) 

00 -b l 2 (k) + ir] 

0 

—^23 



0 oj 

- e3(k) -b irj 

0 



— t23 

0 

00 + e 3 (k) + ir] 

/ 


r 


For T,(k,uj) ~ T,h as in previous case we have 


G(fc,a;)i2 = Si{k,u!)(t23 - -|- e2{k)e3(k)f -b (w - e2(fc))^(w - e3(k))‘^'^ / De{k,uj) 

G{k,00)34 = Si{k,uj) (^{ti 3 t 23 - ti2{Cj - € 3 {k)){ti 3 t 23 + t^iui + e 3 (^))) /De{k,uj) 

G{k,00)53 = Si{k,Lo) {{ti 2 t 23 - ti3{Co - I 2 ik)){ti 2 t 23 + 113(00 + e2(fc))) /De(k,uj) 

De(k,uj) = ti 3 (e 2 (fc)^ - w^) + 4t?2G3^23e3(fc) + ti2{^3(k)'^ - w^) -b 2tl^(tl2( - -b e 2 {k)(ei(k) + Ei(fc, w))) 
-b(w^ - e 2 (fc)^)(w^ + e 3 {k){€i{k) + Ei(fc, w))^ - 4^12^13^23((w^ - t\^)e 2 (k) + uj‘^e 3 (k) + (ei(fc) + Ei(fc, w)) 
x(lo^ — ^23 -b e2{k)e3(k))) + 2t\2{t23i~^^ + i^iik) + Y,i{k,oo))e3(k)) + t\^{—Lo‘^ + 2^23 + ^2(k)e3{k)) -b (lo^ 

-b(ei(fc) -b Ei(fc,a;))e 2 (fc))(w^ - e 3 (k)‘^)) + {(Si(k,u;)^ - w^) -b (ei(fc) -b Ei(fc, w))^)(t 23 - 2 t 23 (w^ 
+e2(k)e3(k)) + (oo^ — e2(k)‘^)(oo'^ — esik)'^))] 


The tunneling strength tij depends on the overlap inte- for G{k,oo). With this approximation, we carry out the 
grals, i.e. the symmetry of the respective eigenstates. 

For the special case of ti 3 = 0 (such as the trilayer case, 
where there is no direct tunneling between the first su¬ 
perconductor layer and third metallic layer) but finite 
and small ti 2 and t 23 , we take the leading order in tij 
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frequency integral and obtain: 


(ct - ) 
' felt -felt' 


(ct - ) 

' fest -fe 2 t' 


(ct c^ - ) 
' feat -fcat' 


-Si{k,U}) 

2Ei{k) 

__ 

2E,{k){E,{k)^-e2{k)^) 


(A18) 

(A19) 


_ -Si{k,uj)tl2tl3 

2E3lk){E3{kY -l2ikY){Ei{kY 



Similar to the bilayer case, we obtain the rough es¬ 
timate for proximity induced gap magnitude on the 

third layer: IA 3 I < | Efe 
.2 .2 

_ ^ 12^23 _I A I 


The same formulation can be extended to the case of 
two bands in a single metallic layer. In this case, we per¬ 
form the same calculations but with ^23 = 0 and nonzero 
but small ti 2 and ^ 13 . The label 2 and 3 here denote the 
two bands in the metallic layer. We will not present the 
calculations here, but just mention that it is a straight¬ 
forward extension of current formula. 


3. Metallic layer with spin orbit interactions 

For two dimensional electron gas in a symmetric quan¬ 
tum well, such as the heterostructure in the bilayer sys¬ 
tem discussed above, or the surface state of three dimen¬ 
sional topological insulator, the anisotropy or asymmetry 
in general leads to some type of spin orbit interactions. 
To include this factor, we consider a generic linear mo¬ 
mentum dependent spin orbit coupling in our two dimen¬ 
sional metallic layer: 

Eso — OlkxCy -\- ^ky(7x - 2^^ “b kyCTx'} 

+ - P){kxCry - kyax) (A21) 

The term times (a + /3) is named Dresslhaus effect, 
often occurring in systems lacking reflection symmetry. 
The other term associated with {a — 0) is called Bychkov- 
Rashba effect which happens when inversion symmetry 
is broken. To explicitly incorporate the spin degree of 
freedom we use four component Nambu bases on a sin¬ 
gle layer: . Here 

I = 1,2 denotes layer index and this basis is chosen as the 
basis for direct product of spin and electron hole space. 
For real Ai, the matrix form of the Hamiltonian describ¬ 
ing the superconducting layer {Hs), normal metal layer 


and tunneling term between the two Elt, in this 
four component basis are: 

Hs = (ei{k) - fii^ Tz (g) CTo + AiTj, g CTo, 

Hn = - fj.2^ g CTo + "Tz g {akxCJy + Pkyax) , 

Ht = tT;,^ao- (A22) 

Note that we use e^\k) to denote the diagonal (in 
spin space) part of the iJjv, and e 2 (k) — ^2 = ^ 2 ^\k) ± 
\J\ctkx\'^ + \l3ky\'^ — fi 2 is used to denote the eigenvalue 
of H]s[ as in other sections. Furthermore, in Ea. (IA22p 
we have assumed Ai = A{k) = A{—k) by placing the 
conventional even pairing superconductor (s or d wave 
superconductor) on top of the metallic layer. The 4x4 
matrix forms of Ea. (|A22|) are: 


/ ei(/) 

0 

Ai(fc) 

0 \ 

0 

ei(fc) 

0 

Ai(-fc) 

Alik) 

0 

-ei(fc) 

0 

V 0 

Ali-k) 

0 

-eiik) J 

^ e2(fc) 

esoik) 

0 

0 \ 

esoik)* 

e2ik) 

0 

0 

0 

0 

-62(4 

-esoik) 

1 0 

0 

-esoik)* 

-e2ik) / 


Here eso{k) = I3ky + iakx- The tunneling term Elt with 
constant tunneling amplitude t takes the form: 


/too 0 \ 

0 t 0 0 

0 0 -t 0 

V 0 0 0 -t j 


This 4x4 matrix for Ht connects the bases of 4'^: and 
'^fc 2 ’ the assumption of momentum conservation. 

In this paper Ai(—fc) = Ai(A:) as our superconductor is 
either s wave or d wave type. The full Hamiltonian, given 

by 


H = 


Hs Ht \ 

Ht Hm ) ’ 


is described by a 8 x 8 matrix or in the (4>j: 4'^ 2 ) bases. 

Assuming t is real and denoting oj = to + ir] as before, 
the matrix element of the retarded Green’s function G = 
{ui — H + related to second layer pairing amplitude 
A 2 are 
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G{k,u})57 =-G{k,u})e8 (A23) 

= Ai{k)t‘^ - 2f + ei(fc)e 2 (^)) + (|Ai(fc)p - + ei(fc)^) (|afca:P + |/3fcyP - + e 2 (^)^)) /Dei{k,uj), 

G(fc,a;)58 = 2A(k)t^{iakx + Pky) {t^ei{k) - (^|Ai(fc)p - e 2 (fc)) /Dei{k,uj), (A24) 

Dei{k,uj) = - 2t^uj^ + (A(fc)^ - A^)(|afca;P + l/S^i/P - + 4i^ei(fc)e2(fc)(G(fc)^ + A(fc)^ - (D^) 

x(|afca;P + |/3fcyP + — e2(fc)^) — 4t^ei(fc)e2(fc)(t"^ — 2t^Cj^) + (A(fc)^ — uj^)e2{k)^{2{t‘^ — 2t^uj^ — (A(fc)^ — w^) 

x(|a/cxP + |,8fcyP + - e2(fc)^))) + ei(fc)^(|aA:a:P + - e2(fc)^)^ + 2ei(fc)^((|afca:p + |/3fcyp)(-t^ 

+A(fc)^(|a/cxP + |/3^yP)) ~ (2A(fc)^(|aA:2;p + |/3A:yp) + {t^ + |afca;p + |/?fcyp)^)w^ + (A(fc)^ + 2(|a/cxP + 
xCj'^ — w® + (3t'^ — 2t^w^ — 2(A(ft)^ — w^)(|afca;p + |/3fcj/P + u)^))e2(fc)^ + (A(fc)^ — w^)e2(fc)'^)- 


This G{k, 10)57 = —G(k,uj)e 8 has to do with the fact that (which is not the case in the main text). Following the 
there is no az term in this spin orbit metallic Hamiltonian same Hartree and small t approximations as in Appendix 

ixn we get 

I 

_ -Ai(fc)t^(|afc.p + \l3ky\^ + e-2(fc)^ - F;i(fc)^) _ ^^25) 

m -k2i 2F;i(fc) ((F;i(fc)2 - \ak,\^ - |/3fc„|2)2 _ 2{E,{k)^ + |afc,|2 + \f3ky\^)e2ik)^ + € 2 (^) 4 ) 

c’f ^ -_ Ai(fc)t^(mfc^ + /3A:^)e2(fc) _ ^^20) 

m -m - ((F;i(fc)2 _ \ak ,\2 _ \pky \^)2 _ 2 {E,{k)^ + \ak,\^ + |/3fc,|2)e2(fc)2 + £ 2 (^) 4 ) 


{ot,P) 

A2s 


Ai 

(1,±1) 

0.07782 

-0.07747i 

0.92121 

-H 

0.09534 

-0.08590i 

0.92104 

(0,0) 

0.14964 

0.0 

0.92181 


TABLE II. Pairing amplitude for various (a, P) with 712 = 
1.5, t = 0.1, Fi = —0.07, and 160 x 160 k-points is used 
in momentum summation. e^\k) = 0 for the top two rows 
and £^\k) = \k\ = + ky for the bottom row. Pairing 

amplitude A 2 p+iy is related to A 2 p+ip listed in the table by 
A2p+ip = A 2 p+ip{asm{kx) + iPsm(ky)). 

From the numerator of Ea. (jA26p and Eo. (jA25p . it is 
clear that the leading pairing amplitude of spin triplet 
(tt or -I'D is of p ± ip form, and that of spin singlet is 
of s or d -wave form (depending on the source super¬ 
conductor, i.e. the symmetry of Ai(fc)). This is indeed 
verified numerically as shown in Table |TT1 and Table IIIII 
in which we use s-wave superconductor as first layer with 
dispersion e“i(/c) = cos{kx) + cos{ky) + 1. In this compu¬ 
tation we have solved the superconducting gap equation 
self consistently, and before coupling the metallic layer 
the gap magnitude of the s wave superconductor is 1 by 
choosing Vi = —0.07. 

Eq. (IA26I) also illustrates why the d-wave pairing in the 
superconductor change the pairing orientations as the dif¬ 


{a,P) 

A2s 

^ 2 p 

Ai 

(1,±1) 

0.02927 

-0.00988* 

0.92093 

-H 

0.03055 

-0.00860i 

0.92085 

(0,0) 

0.30395 

0.0 

0.92009 


TABLE III. Pairing amplitude for various {a, /3) with 
(-^\k) = 1.8 — (cOs(fca;) -I- cos(fcp)), fl 2 = 1.5, t = 0.1, 
Vi = —0.07, and 160 x 160 k-points is used in momentum 
summation. Pairing amplitude A 2 p+ip is related to A 2 p+ip 
listed in the table by A 2 p+ip = A 2 p+ip{cism{kx) + ifism{ky)). 

ferent signs of d wave pairing gap modifies the relative 
sign along k^ and ky directions. Note that there is no p- 
wave pairing in this spin triplet sector (ti + it) but only 
p±ip pairing in the (tt) or (tt). This has to do with the 
lack of (Jz terms in the model Hamiltonian Ea. (jA22p . By 
introducing Zeeman field terms^^ we can also generate 
the p-wave pairing in this (tt + tt) spin sector as shown 
in the discussion of warping terms in Section |TT1 


Appendix B: Comparison between bilayer and 
trilayer system 

Here we compare the proximity induced pairing am¬ 
plitudes of the bottom metallic layer in the bilayer and 
trilayer system. The bottom metallic layers of both sys- 
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terns have the same linear dispersion Eiq = |fc| and 
chemical potential = 1-5- For both bilayer and tri¬ 
layer systems, we choose the same dispersion Esc = 
cos fca;-|-cos ky and chemical potential ^sc = —1 in the su¬ 
perconducting layer. The top and bottom Fermi surfaces 
are chosen to take different forms to highlight the issues 
of mismatched Fermi surfaces, leading to poor supercon¬ 
ducting proximity effect, in general, compared with the 
matched ones. An intriguing question is whether the mid¬ 
dle metallic layer could serve as a good bridging layer to 
enhance this proximity induced pairing amplitude in the 
bottom layer, compared with the direct coupling in the 
bilayer system. 



- 3 - 2-10 1 2 3 




(a)Bilayer FS 


(b) Trilayer FS 


(c)Trilayer FS2 





(d)Pairing amplitude v.s. normalized hopping 


FIG. 4. Ajv = |A:| and Esc = cosfcj, -I- cosky with /tjv = 1-5 
and use = —1.0. Kci = 0.08. 40 x 40 k-points. (a). Fermi 
Surface of the bilayer system (b) .Fermi Surface of the trilayer 
system for Ehi = 4 cos cos ky -\- 2 cos 2kx + 2 cos 2ky with 
fibi = 0.0. (c).Fermi Surface of the trilayer system for Ef,2 = 
(cosfcj, + cosky — 1.05)(cos fca; -I- cos fcji -I- 1.05) with fj,b 2 = 0.0. 
(Green for the SC layer, red for the N layer and blue for the 
bridge layer) (d). Pairing amplitude of the bilayer and trilayer 
systems. 

The tunneling strength of the bilayer and trilayer sys¬ 
tems are labeled by tb and tt, assuming the tunneling 
strengths to be identical in the trilayer system for the 
ease of scaling argument. Compared to the bilayer sys¬ 
tem, electrons in the trilayer system have to hop twice 
from the top superconducting layer to the bottom metal¬ 
lic layer. To compare the results of the trilayer systems 
with those of the bilayer system, we label the effective 
hopping term of the trilayer system as and that of bi¬ 
layer system as tb (Here we choose Ag = 1 so the dimen¬ 
sionless quantity is (tt/Ag)^ = Same for tb/Ao = tb-). 
We then compare the pairing amplitudes of Ni and N 2 
as a function of effective hopping terms tb and t^ respec- 



-3-2-10 1 2 3 



-3 -2-10 1 


(a)Bilayer FS (b)Trilayer FS 



(c)Pairing amplitude v.s. normalized hopping 


FIG. 5. En = |fc| and Esc = cosk^ + cosky with = 1.5 
and use = —1.0. Vsd = 0.08. 40 x 40 k-points. (a). Fermi 
Surface of the bilayer system (b). Fermi Surface of the trilayer 
system for Ebs = |A:| with fibs = 2.2. (Green for the SC layer, 
red for the N layer and blue for the bridge layer) (c). Pairing 
amplitude of the bilayer and trilayer systems. 


lively. 

To test this bridging layer idea, we choose two spe¬ 
cial energy dispersions for the middle metallic layer 
Ebi = 4 cos kx cos ky + 2 cos 2kx + 2 cos 2ky and Eb 2 = 
{coskx+cosky — 1.05){coskx + cosky + 1.05). The chem¬ 
ical potentials fibi = /262 = 0 are chosen such that the 
Fermi surface covering both that of the top and bottom 
layer, as shown in Fig |4(b)| and Fig |4(c)[ We plot the 
pairing amplitude strengths A 2 as a function of the effec¬ 
tive hopping terms tb and t^ for bilayer and trilayer sys¬ 
tems in Fig |4(d)[ As shown in Fig |4(d)[ for small effective 
hopping terms, this choice of bridging layer does enhance 
the pairing amplitude of the bottom metallic layer. For 
larger effective hopping terms, however, the bilayer sys¬ 
tem always has larger pairing amplitudes than the tri¬ 
layer system regardless of the shape of middle metallic 
layer Fermi surface. 

In contrast to the previous case where the Fermi Sur¬ 
face of the bridging layer matches those of the supercon¬ 
ducting and bottom metallic layer, we further study the 
case where the Fermi Surface of the bridging layer lies 
between those of the two layers. The same dispersions 
are chosen for top and bottom layers while the bridging 
layer has a linear dispersion Ebs = |A:| with the chemi¬ 
cal potential fibs > l^N as shown in Fig|5(b)[ The result 
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is shown in Figl3 In contrast to the matching of both 
Fermi surfaces shown in Fig 01 the enhancement here is 
relatively large in terms of magnitude and the range of 
the effective hopping terms. However, at larger tunneling 
strengths the direct coupling in the bilayer system still 
gives larger pairing amplitude. 

All of these numerical results can be understood quali¬ 
tatively from the analytic pairing amplitude derived per- 
turbatively with Hartree approximations in Appendix 
El Under the small tunneling strengths approximation 
the proximity induced superconductivity pairing ampli¬ 
tude for second and third layer in trilayer system are 
given in Ea. (IA19p and Ea. (IA20p . which we rewrite as: 


(ct - ) 

' febt -khV 


(cl - ) 
' fett -ktV 


__ 

2ES){Ei{kY-ei{kY) 

_-^i(fe)4_ 

2ES){Ei{kY - ImikY) {Ei{kY - lb(kY) 


In above, we replace the tunneling amplitudes by 
the bilayer and trilayer ones. The form of the sec¬ 
ond layer pairing amplitude (cl c^ - ) in Ea. dAlQI) is 

k^t —k2i 

the same as the bottom metallic layer pairing ampli¬ 
tude in bilayer system under this Hartree and small 
tunneling strength(s) approximation, and we relabel the 
(cl c^ -, ) by (cl c^ - ). Within these approximations 

k2\ —^ 24 ' 

which is in general smaller than one even if we set A = tb- 
The way to make this ratio larger than one is when 
\Ei{k)\ ~ \em{k)\, the so called resonant condition in 
the scattering theory. This condition is made when the 
middle layer dispersion em{k) are not completely in line 
with the original one given by ei(fc), but with the Bogoli- 
ubov quasi particle level Ei{k). This is why the range 
and magnitude of enhancement shown in Figl5]is greater 
than that in Fig01 At larger tunneling strength, the 
shift of energy level by the single particle tunneling be¬ 
comes more important, and the proximity induced gap 
actually decreases as discussed in AppendixlBj As the 
suppression of pairing amplitude in trilayer is twice of 
the bilayer, the bilayer always has greater pairing ampli¬ 
tude at larger tunneling strengths as shown in Fig0]and 
FigEl 

Knowing that adding an additional middle metallic 
layer in general does not help improve the proximity ef¬ 


fect, we discuss the issues of lattice mismatch, which is 
quite general, but difficult to compute for incommensu¬ 
rate lattice ratios, for interface between different materi¬ 
als in the next section. 


Appendix C: Lattice mismatch issues 


For different materials, the lattice sizes and shapes are 
usually different. Here we focus our discussions on the 
c ases of square lattice, and the lattice misrnatch discusse d 


Cases 

FSM, LM 

FSM, LmM 

FSmM, LM 

FSmM, LmM 

A 2 

-0.42214 

-0.27422 

-0.13747 

-0.25523 

Ai 

0.90649 

0.91302 

0.91354 

0.91299 


TABLE IV. Terminologies used in this table: FSM = Fermi 
surface matched; LM = Lattice matched; FSmM = Fermi 
surface mismatched; LmM = Lattice mismatched. For lattice 
mismatch cases we choose the ratio of lattice size in supercon¬ 
ductor/metal as 3/4. All other model parameters are specified 
in the text of Appendix o 


here means different lattice lengths. We compare the ob¬ 
tained pairing amplitude A 2 in the metallic layer and 
the self consistent gap magnitude Ai of the supercon¬ 
ducting layer. For matched Fermi surfaces and lattice 
sizes we choose ei(fc) = —€ 2 {k) = cos{kx) + cos{ky), = 
—^2 = — I.O, with 1 and 2 denoting superconducting and 
metallic layer and tunneling amplitude t = 0.1. For mis¬ 
matched lattices we choose ei(fc) = cos(3A:a;) -I- cos{3ky), 

Ml = -M 2 = -1.0, -e 2 (fc) = cos(4fca,) -I- cos(4A:y), 
t = 0.1 and Vi = —0.07. Similar choice is done for 
the mismatched Fermi surface with parameters chosen 
as ei{k) = e 2 {k) = cos{kx) + cos(fcy), mi = —M 2 = —1.0, 
t = 0.1, and Vi = —0.07. The results are listed in Table 

EyI 

It is clear that for originally matched Fermi surfaces, 
the factor of mismatched lattice makes less overlapping 
region of Fermi levels. Therefore we shall expect the 
decrease in the proximity effect. For the originally mis¬ 
matched Fermi surfaces, this lattice mismatch factor ac¬ 
tually increases the percentage of overlapping Fermi lev¬ 
els, and possibly enhance the proximity effect compared 
with original identical lattices. These intuitions are con¬ 
sistent with what we obtained in the numerical results 
shown in Table EYI 
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